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Abstract 

Gravitational waves propagate along null geodesies like light rays in the geometrical optics ap- 
proximation, and they may have a chance to suffer from gravitational lensing by intervening objects, 
as is the case for electromagnetic waves. Long wavelength of gravitational waves and compact- 
ness of possible sources may enable us to extract information in the interference among the lensed 
images. We point out that the interference term contains information of relative transverse ve- 
locity of the source-lens-observer system, which may be obtained by possible future space-borne 
gravitational wave detectors such as BBO/DECIGO. 
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I. INTRODUCTION 



Direct observations of gravitational waves from astrophysical sources would open a new 
window to the universe, and a new era of astronomy. It is now well known that there are 
many gravitational lens systems of various sorts in the electromagnetic astronomy. This 
must also happens in gravitational wave astronomy in the future. One difference between 
electromagnetic waves and gravitational waves are their wavelength. The wavelength of 
gravitational waves which the currently working and the planned gravitational wave obser- 
vatories may detect are typically much larger than the electromagnetic waves in the current 
astronomy. Therefore once they suffer from gravitational lensing, we have a chance to ob- 
serve interference. To be more quantitatively, a condition may be met for interference to be 
seen. In the electromagnetic astronomy, a wave source is typically extended and it consists 
of pieces of incoherent emitters. Then, to observe interference, the typical linear dimension 
of the source / should be smaller than or comparable to the wavelength A divided by the 
angular separation among the multiple images 9 (in radians), or I ^ k,X/6 = kIq with n of 
order unity in a typical gravitational lensing occurring in cosmology. Otherwise incoherent 
waves from different parts of the source distract the interference. Observing at, say, A = lyum 
and assuming 6 = 10" ~ 5 x 10^^ radians, Iq — 2 cm. This is one reason why we usually 
do not see interference in the electromagnetic astronomy. One exception where wave effects 
become important is scintillation for radio sources due to ionized interstellar /interplanetary 
gas (See, e.g.. Chap. 13.4 of [Ij). As an application, for example, the size of a gamma- 
ray burst (GRB) can be estimated using the interstellar scintillation of radio afterglow of 
a GRB (See, e.g., On the other hand, there are compact gravitational wave sources 

whose linear sizes are of order the wavelength they emit. These include relativistic coalesc- 
ing binaries such as neutron stars binaries and isolated non-axisymmetric pulsars, both of 
which are two of the most promising sources (e.g., j^]) for the gravitational wave detectors 



that are currently working or planned, such as LIGO '4'] , LCGT [5] , LISA [6^, BBO (7|, and 



[el, BBO fzi, 

DECIGO [8]. At around the BBO/DECIGO observing frequency of / = O.lHz, the orbital 
separation / of a neutron stars binary would be I ~ lO^km(2MNs/2.8M0)^/'^(//O.lIIz)~^/^, 
while /o ~ 4 X 10'^A[J{f/0.mz)-\e/10")-^ » I. 
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In the above argument, we saw that there is an upper hmit on the size of the emitter 
to see interference. This is not a sufficient condition for interference pattern to form in the 
space and for us to observe its amphtude changing as the Earth (or a detector) moves in that 
spatial pattern. For a clear interference pattern to form in the space, coherence time must 
be larger than the time delay among multiple images. Note that the coherence time is of 
order 1/A/ where A/ is the frequency shift during an observation. Suppose an extragalactic 
neutron stars binary as a gravitational wave source. Then to have a large enough signal to 
noise ratio with BBO/DECIGO to claim detection, we may integrate the signal from the 
binary for a couple of years during which the signal frequency of the source binary would shift 
over of order of the BBO/DECIGO observing frequency, / ~ A/ ~ 0.1 Hz. The coherence 
time would then be 10 seconds, which is smaller than, for example, the time delays of the 
known gravitational lensing systems listed in the CASTLE data base [9|, HOj. So unless the 
time delay is smaller than ~ 10 seconds or the source is highly monochromatic, we are not 
in a situation where we (on the Earth) move in a clear spatial interference pattern. 

The above argument does not preclude a possibility of detecting an interference term in 
the time domain, or equivalently, in the frequency domain. In this paper, we do not aim to 
detect an interference pattern in the space. Rather, we propose a method of a filtered cross- 
correlation of gravitational lens images in the frequency domain. We shall point out that the 
interference term contains information of the relative transverse velocity of the source-lens- 
observer system and study how well we could determine it using the future BBO/DECIGO 
detector. 

There have already been several works that point out importance of interference and, 
more broadly, wave optics effects in gravitational lensing phenomena of gravitational waves 
14 . fisj l . Since a gravitational wave source such as a coalescing binary is essentially 



liiU. 

Q 
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131, 



a point source, we have to use wave optics rather than geometrical optics approximation for 



such a source near caustics 
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17|. The works [ij, |l5| studied observable effects on the 



gravitational waves propagating in an inhomogeneous universe. In this paper, we propose a 
method for extracting an interference term in the geometrical optics limit and for obtaining 
information on a relative transverse velocity of the source-lens-observer system. There are 
various methods to measure a relative transverse velocity of astronomical objects at a cos- 
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mological distance in the electromagnetic astronomy. Those includes for example, a method 
to measure transverse velocities of clusters of galaxies detecting a signature of gravitational 
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and to measure those of galaxies 



21 



22]. Possibilities are also dis- 



scattering of the cosmic microwave background [18|, 
using a parallax effect due to gravitational micro-lensing 
cussed in in the gravitational lens literature by 23|]. Gravitational lensing of gravitational 
waves possibly offers another way to measure a relative transverse velocity of cosmological 
objects in addition to those methods above. 

Our method could be applied to various sources including continuous waves from pulsars, 
cosmic strings and so on. However, to be specific, we consider a coalescence neutron stars 
binary at a cosmological distance. 

This paper is organized as follows. In the next section, we explain our situation by 
discussing a (astronomical scale) Young's interference experiment. In Sec. Ill, we remind 
readers of the basic equations and the notations for study of gravitational lensing of gravi- 
tational waves, for which we follow {2^. After remarks on a situation we would be in when 
searching for gravitational lensing of gravitational waves and on the geometrical optics ap- 
proximation. Sec IV explains how to extract interference term when the approximation 
applies. In the section |Vl we propose a method for extracting information contained in the 
interference term, namely, the relative transverse velocity of the source-lens-observer system. 
We then study a correlation between the time delay and the transverse velocity in Sec IVII 
and a correlation among parameters in IVIII in the absence of noise. Sec. IVIIII then shows 
the precision to which we could measure a transverse velocity by introducing detector noise. 
The section HXl gives summary of our result. 

In this paper we use a unit of G = c = 1 unless otherwise written explicitly. An amplitude 
of a gravitational wave inversely depends on the luminosity distance from the observer to 
the source, to compute which we assume a flat universe with = 0.7, i^matter = 0.3, and 
h = 0.7. 
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II. YOUNG'S INTERFERENCE EXPERIMENT 



To explain our situation, consider the (astronomical scale) Young's two pin-holes inter- 
ference experiment. The (gravitational) wave intensity / at the observer's time and position 
{t, x) on the screen is 



I{t,x) 



ri r2 



' 1 

~2-(l + COs((7i-(72)), (1) 



where is the distance from the k-ih. pin-hole to the observer and is the phase of the 
wave from the /c-th pin-hole. We neglect the difference in ri and r2 in the amplitude and 
write ri ~ r2 = r. The phase difference ait^x) = Ui — may be written as 

a{t,x) = 27r f{t')dt', (2) 

Jtr 

where tr — tr{t, x) — t — Ti{t, x) is the source retarded time for the wave through the pin-hole 
1, Tk{t,x) is the time length for a wave crest to propagate from the source to the observer 
through the k-th pin-hole, and Atd{t,x) = T2{t,x) — Ti{t,x) is the time delay between two 
pin-holes. 

For a monochromatic wave at the frequency f{t) — fc — constant, a{t, x) — 277 fcAtd{t, x) 

and the scale of the interference pattern is Ax — c/9/ f/2 {9 being the angular separation 

of the two pin-holes seen from the observer) as usual. Known possible monochromatic 

gravitational wave sources are isolated pulsars. The planned next generation gravitational 

wave antennas would see those within the Milky way galaxy. So considering micro-lensing 

phenomena on a gravitational wave emitted by an isolated millisecond pulsar in a globular 

cluster (say), the spatial scale of the interference pattern is 

in / lmas\ / lkHz\ , , , 

Aa; ~ 3 X lO^^km ( — ^ ) \ J " 200AU. (3) 

Suppose, in the spatial interference pattern due to the monochromatic wave, we move at 
200 km/s (about the orbital speed of the Solar system around the Galactic center), it takes 
~ 5 years for us to go across the spatial interference pattern. The wave intensity will vary 
in time as 

I{t, x) ~ cos ^TT^^t + const. ^ -I- const. (4) 



In principle, for a monochromatic wave, we would be able to detect relative transverse 
velocity of the source-lens-observer system by using interference pattern (by measuring the 
temporal variation of I{t,x)). If we recognise such interference pattern, we would conclude 
that gravitational force shows interference phenomena and propagate as "wave" . Other than 
that, the information we could obtain, v6, is in principle the same as that in micro- lensing 
experiments in the conventional electromagnetic astronomy. (And the number of millisecond 
pulsars is much smaller than that of the stars in the Magellanic clouds, say.) 

When observing pulsars within the Milky way, it is also important to note that diffraction 
effect becomes non- negligible when gravitational wavelength A is larger than the lens mass 



Mr 



2J|, or Ml ^ lO^M0(//lkHz) ^. If we observe interference of gravitational waves from 



pulsars in the Milky way, then the lens should have mass larger than ~ lO^M© 

Eq. ([3]) suggests a cosmological application of gravitational lensing phenomena on grav- 
itational waves. 

V— (^)(T)(^)- 

Unfortunately for our purpose of seeing an interference pattern, the most promising cosmo- 
logical sources, relativistic compact binaries in an inspiralling phase, do not emit monochro- 
matic wave. For such a source the frequency of the gravitational wave varies in time, to the 
lowest order, 

dt 5 \ J ^ ' 

with Mc = (mim2)^^^/(mi + 1112)^^^ is the chirp mass, being the mass of the k-th. star 
in the binary. For wave having time-varying frequency df{t)/dt = Af^{t) with A and rj 
constants, 



a{t, x) 

271 



iiv - l)(tc - t - At,(t, x))A + f^-J)^ - iiv - l)(te - t)A + f^-J)^, (7) 



with f{t = tc) = /max at which inspiralling phase ends and the binary stars crush into each 
other. The interference pattern changes with time even if there is no relative transverse 
velocity in the system. Moreover because of the smallness of the amplitudes of possible 
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gravitational waves, we normally need to integrate signals for some time duration to get a 
sufficiently-large-for-claiming-detection signal to noise ratio. Schematically, we would do 



Tobs is determined by either the source lifetime, its pass-time over the frequency band of 
our detector, or a pre-set threshold for claiming detection (, say, signal to noise ratio larger 
than 5). Large Tobs in the integration makes the interference term vanishingly smaller than 
the bolometric (non-interference) term. For a massive compact binary, depending on binary 
masses, it is possible that Tobs << (time delay) so that one wave packet through the pin- 
hole 1 comes to the detector, and after some quiet period, the second through the pin-hole 
2 visits it. In this case, no interference pattern forms in the space and we would detect 
their interference by taking a cross-correlation (in a computer). This paper shall use a 
filtered cross-correlation technique to take care of the issues above. But before moving onto 
explaining the technique we use in this paper, we start our discussion by explaining lensed 
waveform of gravitational waves in the next section. 

III. LENSED WAVEFORM OF GRAVITATIONAL WAVE 

We consider a coalescence compact stars binary at the redshift zs with the redshifted 
chirp mass M^.^ = Mc(l + zs) as a target source to be lensed. The detector is assumed to 
be a space borne interferometer such as (one-triangle) BBO/DECIGO, which outputs two 
independent data streams ha=iji{f)- When the gravitational wave from the binary is lensed 
by a lensing object of redshifted mass M^^ = Ml{1 + z^), the wave through the j-th image 
to the detector in the frequency domain would be written, in the geometrical optics limit 
and in the stationary phase approximation, as 24 1 




(8) 




(9) 



with 



$a,j(^) = 0D,i(t) +0pj>(t), 



(10) 
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where nj = 0,1/2 and 1 when the j-th image corresponds to a minimum, saddle, and 
maximum point, respectively. Aa,j{t) is written in terms of the detector beam pattern 
functions F^{t),F^{t) [25] and defined as 

A.,,(t) = {{l + {L.N,ffF^{t)+A{L.N,fF:{t)y\ (11) 

where L (given by Ol,4>l) and Nj (given by 6s,j,(t>s,j) are the direction unit vector of the 
binary orbital angular momentum and the direction unit vector toward the j-th image. These 
vectors are defined in the fixed barycenter frame of the solar system. The detector phase 
4>D,jit) and the source phase (ppj^ait) are 

(f)D,jit) = 27r/(t)i?sin% cos (0(t) - ^s,j) , (12) 
0.,.(t) ^ tan-^ f 2(r .i^)F.x(t) \ 

where 0(t) = 2nt/T + ^o with R = lAU and T = 1 year for BBO/DECIGO (Those detectors 
are planned to orbit around the solar system barycenter with its average orbital radius 1 
AU. For simplicity we assumed a circular orbit for the detector here). The amplitude of 
the gravitational wave A and the phase of the gravitational wave ^E'j(/), to the Newtonian 
approximation, are 

^,if) = 2nfU, + ^if) = 2nft,, _ 0^ _ ^ + ^ {SnM^jy'/' , (14) 
1 / 5 ^ 



A = —-^ - 7r-2/3M,5/6 (15) 

where Ds is the source angular diameter distance, fij and t^j are the amplification factor 
and the time delay for the j-th image measured with respect to the fictitious time of arrival 
of non-lensed signal, tcj = tc + tdj- tc and (pc are the time and the phase of the binary 
coalescence. Finally, the time variable t in the above equations must be understood to be a 
function of frequency as 

t{f) = = - 5i87rf)-'/'M-'/', (16) 

W) = t{f) + = te, - 5{8nf)-'/'M-'/\ (17) 

In a gravitational wave search, the unknown parameters would be, to the lowest order, 

{Mcz,(pc,tc,Ds,9s,j,4>S,j,L,hil - K2, Atd = td,2 -td,l, |/il|/|/i2|}, (18) 
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where the last three depend on the lens property. 

IV. SEARCH FOR A LENSED GRAVITATIONAL WAVE 

We consider the following situation. First of all, for simplicity we study that two images 
occur due to gravitational lensing. Then the wave from one lens image would be in the 
detector's observation frequency band in the time duration from, say, ts^i to te,i- The 
second wave would be in the band after the time delay At^ later from 2 = + to 
te,2 = te,i + Atrf. We assume the observation frequency and the time delay satisfy fAtd ~ 
0.1(//0.1Hz)Atrf >> 1 so that we can use the geometrical optics approximation ([24] studied 
the case where wave optics is important). As as result, we may detect the two waves 
separately using unlensed waveform templates as if those were unlensed, in the same sky 
direction (The angular resolution of the gravitational wave source is ~ 10~^(Signal to noise 
ratio) ~^ radians with the current BBO design \2(y!\, so we do not expect to resolve images). 
We would then realize that those two waves have similar parameters sets Mcz, (pc L, and 
conclude that those two are actually due to one lensed event. We may find the (relative) 
time delay Atd = — td,i and the ratio of the magnification |/Ui|^''^/|/i2|^''^ = ^iM2- The 
observational errors of these two numbers are at the same level as those of tc and A. We might 
now cross-correlate the two waves to see interference between these two (in a computer). 
However, a simple cross-correlation disappears when the geometrical optics approximation 
applies (and this is the reason we could detect two waves separately, anyway). Let us see 
this using some equations below and propose a method to extract information contained in 
the interference term even in this case. 

A. How to extract the interference term when the geometrical optics approxima- 
tion appHes 

We consider a sort of cross-correlation between the two waves, given by iQ. We use the 
waves with the supports in the time domain assumed to be from te,i — At^ to te,i for the 
wave that reaches the detector first, and from te,2 — Atd = te,i to te,2 for the second. The 
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detector's output Sk(t) in the two time segments {k = A, B) are 

^Ait) = ^i(^) ^'^'l - Atrf < t < te,u (19) 

j 

S^(t) = h2{t) for < t < te,2- (20) 

The Fourier expansions of the gravitational wave hi{t) in the time segment A and that of 
h2{t) in the time segment B have a support from to fe in the frequency domain. The 
frequency domain support of h2{t) in the time segment A is lower than /j . We now compute 
a simple cross-correlation between the two waves as 

with Sh the one-sided spectral density of the detector's noise. Because we defined the 
integral region from /j to fe so that there is no auto-correlation term of h^2if) which 
otherwise appears due to h2(t) term in Syi(t). Now, h^j{f) is given by Eq. ([9]) and the 
simple cross-correlation becomes 



1/2 1 



7/3 



^ / 4w7T^",i(^i)^"(^2)e'"'^^*'^-'"^"^-"^)e'*'^-2(*^)-**"-^(*i)(i/. (22) 

This integrand is oscillatory and the result is effectively zero if |/At(i| >> 1. 

Now, how can we extract non-zero cross-correlation? We propose a filtered cross- 
correlation statistic. We multiply the two images, and before the frequency domain in- 
tegration we further multiply cos(27r/AtJ + P"^)) where Atd = td,2 — the tuplet 

P = (tc,l,tc,2, Mcz, - K2, A^, L), 

e(/,p) ^ $a,2(t2(/)) - $.,l(tl(/)) + TviK^ - K2), (23) 

and denoting a template. We define our detection statistic C« as 



Ca = 4i?e J s,if) '°"^27r/Ati + e(/, p^ ))df 

rfe f-7/3 

X ^A„,i(ti)A„,2(t2)cos(27r/(At,-AtJ) + e(/,p)-e(/,p^))d/. (24) 
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Maximizing an absolute value of the detection statistic \(\ (where ( = XlaCa); we obtain 
our estimate of the time delay At^i. In the next section we study an effect on our statistic 
of a relative transverse velocity of the source-lens-observer system. In the later sections, 
we shall study precision to which we could measure At^ and the transverse velocity and 
correlations among errors of the parameters At^, the transverse velocity, and p. 



V. TRANSVERSE VELOCITY 

We can measure {t^j, Mcz,0sj4'Sj L,nj} using gravitational lensing phenomena in the 
usual electromagnetic astronomy. When one can use information in an interference term, 
we could get new information on the lens/source object: its relative transverse velocity. 
When there is a relative velocity among the source-lens-observer, there may be two effects. 
One is the time variation of the direction of the images Nj, which may not be detected for 
cosmological sources in the foreseeable future because of the poor pointing ability of the 
planned gravitational wave detectors. The other is the Doppler effect which causes rescaling 
of the mass of the source, the time variables and the amplitudes of the gravitational waves. 
The rescaling due to the Doppler effect differs for different images as 

iAf) = kjk - 5{8nf)-'/'M-f, (25) 

and Mcz,j = 'jjMcz, icjk = Iktcj, and Aj = 7^- "^^^A with 7^ = (1 + Nj ■ 0), where the velocity 
(3 = v/c is a linear combination of the observer velocity vq, the lens velocity vl, and the 



source velocity vs as [27 



28| 



^ = + J—-J^^s J, VL, (26) 

with Dls being the angular diameter distance between the lens and the source. 
The lensed gravitational waveform from a coalescing binary now becomes 



with 



TT 3 / - \ -5/3 

^jkif) = 27rftc,jk - 0c - - + ^ [S7iM,z,kf) . (28) 
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If there is only one image, then we could not resolve the degeneracy and would get, in 
principle, a biased estimate of the masses and the source distance. However, if we know in 
advance that we see multi- waves due to gravitational lensing (see Sec. lIVp . then we could 
find the differences of the Doppler factors among different images. 

For later convenience, we introduce a notation for the ratio of the Doppler factors as 

r = ^ = 1 + (iVi - iVa) ■ /3 + 0{(3^). (29) 

72 

Note that F — 1 ~ 6(3± where 6 is the angular size of the lens. Finding F using the interference 
term of the lensed gravitational wave then gives us the transverse velocity f3± [sil , given that 
6 is known from, say, electromagnetic observation of the lens object (say, galaxies/clusters) 
and the host object (say, a galaxy) of the gravitational wave source. 

We compute a cross-correlation between the lensed wave with itself but with frequency 
shifted by F"^ here denoting a template, and further multiply the cosine filter as in the 
previous section. 



C = 4i?e ^ _^/,^^^(/)/,^*^ (r^/) cos(27r/AtJ,, + e(/, p^))) 

f^' df 



3\^^,nf^,\'/'A,A,Re -^/-^/3Hcos(27r/AtJ,, + e(/,p^))). (30) 



The function 7i is defined as 



n = Aa,l(tll(/))A«,2(t22(r^/)) 

X gi-I'a.2(t22(r^/))-i<I>a,l(*ll(/))ei*ll(/)~i'I'22(r^/)^ ^^l^ 



with 



3 / ~ \ -5/3 

^ii(/)-^22(F^/) = 27r/t,,n + - (87rM,,,i/j 

-27r/F^t,,22 - I (87rF^M,,,2/) • (32) 

The parameters to be searched for are, in addition to the parameters p"^, the time delay 
between images, At^, and the F parameter. When F-^ ~ 1 + 9P± + S{9p±) with 6^/5 j_ << 1 
and 6{ef3±) « 1, 

^ii(/) - ^22(r^/) ^ 27r/At,,i - 27r/t22(/)5(^/5±), (33) 
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with Aid,i = 7iAtd. 

Eq. (133|) tells us that larger the observation frequency / is, and/or longer the observation 
time is, with a better precision we measure the time delay At a and/or F 321]. This latter 
tendency will be studied further in the sections IVIII and IVIIII Eq. (!33|) also indicates a 
correlation between At^i and F, which we will study in the next section. 



VI. CORRELATION BETWEEN THE TIME DELAY AND THE TRANSVERSE 
VELOCITY 

When one uses usual matched filtering technique to unlensed gravitational waves, the 
search parameters are 8 in the lowest order (See the list ( |T8|) ). To find the transverse 
velocity of the source, in addition to those, we search for two parameters: the time delay 
Atrf and the transverse velocity F. Since these two parameters are the new aspects in our 
paper, in this section we show how our detection statistic ( (where ( = J2a Ca) depends on 
these two parameters. 

To demonstrate a correlation between Atd and F, we consider a source with the param- 
eters listed in the table [B The source redshift zs, the lens redshift zl, the image separation 
6, the time-delay and the direction of the source are taken from the famous lens system 
Q0957+561 0, Q]. With those parameters, we compute our statistic ( in the F^-AtJ plane 
with other parameters fixed to be the right values. Fig. [1] shows the contour map of (, which 
clearly shows a strong correlation between F"^ and At J. As can be inferred from Eq. fl33l) . 
the correlation is found well approximated by 

AtJ-Atrf = -(F^-F)Tobs, (34) 

near AtJ = At^^. In fact, along this line and projected onto F, the detection statistic behaves 
as in Fig. [2], which shows the maximum appears at F^ ~ F: The maximum of ( occurs 
slightly shifted F^ value from the true value (~ 1%). This happens because in Eq. ( |30l) . we 
should have used instead of 6(/, p) defined in Eq. (!23l) . 6(/, p) of the form 

e(/,p) = <l'„,2(t2(r/)) - ^aAhif)) + 7i{k, - K^). (35) 
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TABLE I: The source, lens, detector parameters for the figures [T] and [2j Here Tobs is the integration 
time (time for the wave to pass the frequency band from /j Hz to fe Hz) in years, and v = \v\: 
magnitude of the relative transverse velocity of the source- lens-observer system defined in Eq. (|26p . 
Other parameters are explained in the main body of the text. The source redshift zs, the lens 
redshift zl, the image separation 9, the time-delay and the direction of the source are taken from 



the lens system Q0957+561 [9|, 



29|. 



fi [Hz] 


fe [Hz] 


Tobs [yrs] 


Atd [yrs] 


V [km/s] e ["] 


r- 1 


0.113 


1 


1.137 


1.14 


480 6.26 


4.9 X 10^*^ 




tc [yrs] 


mi, 1712 [Mq] 


9l [rad] 


(f)L [rad] 9s [rad] 


4>S [rad] 


1.5,1 


1.14 


1.4,1.4 


1.09 


2.90 2.62 


0.99 


4>o [rad] 


4>c [rad] 


aQ [rad] 


zs 


ZL 




2.64 


2.55 


2.06 


1.41 


0.36 





Unfortunately, since use of this demanded large computational power in our code, we used 
6 instead paying ~ 1% bias in F — 1. 



VII. CORRELATION AMONG THE PARAMETERS 

As is mentioned in Sec. llVt we specifically assume the following situation: Two images 
occur due to gravitational lensing. A gravitational wave is in the detector's observation fre- 
quency band during the time segment from ts^i to te,i- We use matched filtering technique to 
estimate the source parameters (8 in the lowest order, Eq. ( |T8l) ). Another wave reaches the 
detector from ts^2 = ts,i + to te,2 = te,i + '-^td, and we use matched filtering technique to 
estimate the source parameters. From similarity of the estimated values of the parameters, 
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FIG. 1: ( contour map in the — T and At J — At^ plane, in the absence of noise. The source, 
lens, and detector parameters are listed in the table HI All the template parameters other than 
and At J are fixed to be the right values. The dotted line is Eq. (j34p . which approximately follows 
the ridge of C = C(r^,AtJ). 
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we would conclude that we see gravitational lensing of gravitational waves. We then con- 
struct a filter 0(/, p^) from the estimated parameters, compute a filtered cross-correlation 
(, and estimate F and At^. Therefore, in practice, an estimates of F and Atd should be 
affected by errors in estimates of other parameters p"^ in 0(/, p"^) (We do not estimate p"^ 
at the same time as F"^ and At^^ when computing the filtered cross-correlation). 

To find the effect of the errors in p-^ to the estimate of F, we performed the following sim- 
ulations. We generate a wave with the source parameters 9l, (pL, Os, 4>s, 4>c, 4>o, cto randomly 
chosen (and thus we generate waves with randomly chosen p). Other source parameters are 
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FIG. 2: along the line Eq. (|34p projected onto the — T axis, in the absence of noise. The 
source, lens, and detector parameters are listed in the table HI All the template parameters other 
than and At J are fixed to be the right values. 
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fixed and listed in the table [B For each parameters set, we take the template parameters by 
adding random Gaussian errors 6p to the injected signal parameters p as = p + 6p. The 
standard deviations of the errors are 10~^ for the relative errors in the mass rrij and the time 



of coalescence tc, and 0.1 radians for the angular parameters 9l,(Pl,0s,4>s,4>c,4>o^C(o [30 |. 
We then compute our statistic ( and find F and Atd that maximize \(\. We repeated the 
above steps 1100 times. Fig. [3] shows the result of the simulation. Here we search over the 
region of |F^ — F| < lO""^!"! \AtJ — Atd\ < 20 sec (The size of the search region is limited from 
our computer power). We see the estimates of F and Atd roughly satisfy Eq. (!34|) . Several 
points (8% of all) are accumulated around the upper left and the lower right boundaries of 
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the search region, for which the estimates should in fact be outside of the region. Fig. H] 
shows the result of the same simulation but with the true time delay Atd enlarged to 5.1 
years, and correspondingly = 0.065 Hz. We see a better correlation to Eq. and a 
smaller points accumulated around the boundary of the region (3%). Both of the figures 
show a tendency that larger the value of ( is, better the accuracy of our estimates for F and 
Atd are. 

Now what are the effect of the errors in the onto estimates of F and At^? Fig. 
|5] and Fig. [6] show the cumulative probabilities of the absolute errors in |F^ — F| and 
|AtJ — Atd\. From FigJSl we conclude that errors in cause errors of less than 2.5 x 10~® 
in F"^ (F = 4.9 x 10^^) for 35% of the simulation for Atd = 1-14 years and 70% of the 
simulation for Atd = 5.1 years. 

Finally, the scatter plot Fig. [7] shows how the errors in Mcz and tc affect our estimate of 
F. Since Mcz,i/Mcz,2 = we can determine F using estimates of M^zj, if this ratio can be 
determined accurately enough. This figure shows that even though we admit errors in Mcz 
of order 10~^ (note that the injected signals have F = 4.9 x 10^^), F can be determined at 
less than 10^^ accuracy in our simulation, using our filtered cross-correlation method. 



VIII. PRECISION TO WHICH WE COULD MEASURE 6I/3_l 

In reahty, there is a noise in detector's outputs, so the Eqs. fl2Ul) and fl?I]) should be 
actually 




(36) 



i=i,2 
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FIG. 3: A scatter plot of maximum locations of ( in the — T and AtJ — Atd plane for 1100 
sets of randomly chosen source, lens, and detector parameters. Here At^ = 1.14 years. Symbols 
indicate C values: Cmax ~ 2.8 x 10^ denoting the maximum value of ( in the 1100 trials, diamonds 
for < C < 0.25Cmax, crosses for 0.25Cmax < C < O.SCmax, circles for O.SCmax < C < 0.75Cmax, and 



squares for 0.75Cmax < C ^ Cmax- See Sec. IVIII for details. 




where nk=A,B(t) are the detector's noise. The presence of noise affects firstly the estimates 
of the filter parameters when one use matched filtering to estimate those. Secondly, it 
affects the estimates of and AtJ when we use a filtered cross-correlation technique. We 
studied the former effect in the previous section, and here we study the combination of them. 
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FIG. 4: The same as Fig. [3]but here Atd = 5.1 years and correspondingly /j = 0.065 Hz. See Sec. 
IVnifor details. 
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In the presence of the noise, ( statistic should be 

C = 4i?e /'^ '^^P'^y^ cos(27r/AtT,, + e{f,p^M 

= Co + ^Rej: r "^^o cos(2vr/AtI, + 9(/,p^))^/ 



a f'' 



Shif) 



Shif) 



cos(27r/At^_i + e(/,p^))rf/ 



+ ^Ref^^ cos(2vr/At^, + 9(/,p^))c^/. 



(39) 



Co is C statistic without noise given by Eq. fl21|) . As is in the previous section, we first generate 
a wave with a randomly chosen source parameters set l, 4>l,^s, 4>s, 4>c, 4>o,<^o where other 
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FIG. 5: Cumulative probability distributions of |r-^ — r|. See Sec. IVIII for details. 
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parameters are taken from the table [H We then assume template parameters by adding 
Gaussian errors to the true parameters of the generated signal. The standard deviations of 
the errors are the same as before. For each realization of the wave, we further add random 
Gaussian noise nj{f) satisfying < n{f)n*{f') >= Sh{f)S{f — f')/2 where < ... > denotes 
ensemble average. We then compute our statistic ( and find F and Atd that maximize \(\. 
For each generation of a wave, noise frequency series is generated 1000 times. We went 
through these steps 10 times, so that we obtained total of 10 x 1000 estimates of F and Atd- 
The results of the simulation are as follows. The cumulative probabilities of the errors 
in |F^ — Mean(F-^)| are shown for 8 out of 10 waves in Fig. [8] (For the excluded 2 waves 
the maximum of the Co appears to close to the search region boundary and we could not 
compute the cumulative probabilities properly). Here Mean(F^) for each source parameters 
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FIG. 6: Cumulative probability distributions of |AtJ — Atd\. See Sec. IVIII for details. 
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set is the mean of over 1000 realizations of the noise series and does not necessarily equals 
the true value of F, as studied in the previous section. The standard deviation of the errors 
is 1.9 X 10~^ and 90% of the time the error is less than 3 x 10~^ on average of 8 waves. We 
also performed the same analysis but with Atd = 1.14 years replaced by At^ = 5.1 years 
(so that fi = 0.113 Hz by = 0.065 Hz). Fig. shows the result for this longer time delay 
(and thus longer integration time of 5 years). The standard deviation of the errors in this 
case is 1.1 X 10~^ and 90% of the time the error is less than 2 x 10^^ on average of 10 waves. 

Finally, we performed the same analysis for v = lOOOkm/s or F = 10^^, and obtained 
quantitatively the same results as above. Combining the results in the previous sections, 
we conclude that F^ would be determined, in our example, with the future (one-triangle) 
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FIG. 7: A scatter plot of maximum locations of ( in the MJ^/Mcz — ^ and tj /tc—1 plane for 1100 sets 
of randomly chosen source, lens, and detector parameters. Here At^ = 1.14 years. Symbols indicate 
|r^-r| values: dots for < |r^-r| < 3.33xl0-^ crosses for 3.33 x 10"^ < |r^-r| < 6.67x10"^, 
and squares for 6.67 x 10^^ < |r-^ — r|. The true value of T is 4.9 x 10^^. 
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(40) 



(41) 



for 5.1 years integration, where the first error in each equation is due to the errors in the 
estimates of the filter parameters p and the second error comes from the noise terms in the 
filtered cross-correlation analysis in Eq. (l39l) . 
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FIG. 8: Cumulative probability distributions of IF"^ — Mean(r^)| for 8 waves with different source- 
lens-detector parameters sets in the presence of noise. Here At^ = 1.41 years. See Sec. IVIIII for 
details. 




IX. SUMMARY 

In this paper, we considered a gravitational lensing phenomena of gravitational waves in 
the case when the geometrical optics approximation applies. To be specific, the source was 
assumed to be a coalescing neutron stars binary at a cosmological distance and a detector 
to be the planned space-borne detector (one-triangle) BBO/DECIGO. We then proposed a 
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FIG. 9: Cumulative probability distributions of |r^ — Mean(r"^)| for 10 waves with different 
source-lens-detector parameters sets in the presence of noise. Here At^^ = 5.1 years. See Sec. IVIIII 
for details. 
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filtered cross-correlation method of extracting the relative transverse velocity of the source- 
lens-observer system using an interference term. 

We performed series of simulations to study correlations among parameters and to es- 
timate errors due to detector's noise. In these simulations, we adopted lens parameters 
taken from the lens system Q0957-I-561 [ol, [2^ as a reference for a demonstration purpose, 
in which the time delay is Atd = 1.41 years. With our method, we found the absolute error 
in r = 6l3± {6 is the images separation and l3± is the relative transverse velocity), with 
the future (one-triangle) BBO/DECIGO detectors to be 6(T) = ±5.5 x 10"* ± 1.9 x 10"^ 
for the time delay of 1.14 years, where the first error in the equation is due to the error in 
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the estimates of the filter parameters p and the second error comes from the noise terms 
in the fihered cross-correlation analysis in Eq. fl39|) . The errors in the equation are the 
standard deviation estimated in our simulations. We also performed the same analysis but 
with Atd = 5.1 years. In this case, we found 6iT) = ±2.0 x 10"^ ± 1.1 x 10"*^. 

Although the probabilities of lens phenomena is highly uncertain mainly due to the fact 
that the event rate of neutron stars binaries coalescence is highly uncertain, one observation 
of transverse velocity of a neutron stars binary may give us valuable information of such 
population. 
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